clear all
close all
clc


D=load('TimeHistory');

St=2;
ed=length(D);
% ed=5006;
Strain0=18;
Stress0=13;

subplot(1,2,1);
plot(-(D(St:ed,Stress0+1)+D(St:ed,Stress0+2)+D(St:ed,Stress0+3))/3000,(D(St:ed,Stress0+1)-D(St:ed,Stress0+2))/1000,'color','b','linewidth',1.5);

xlabel('p (kPa)','Fontname','Times New Roman','Fontsize',16);
ylabel('\sigma_1-\sigma_3 (kPa)','Fontname','Times New Roman','Fontsize',16);
set(gca,'Fontname','Times New Roman','FontSize',16,'linewidth',1.0);


subplot(1,2,2);
plot(-D(St:ed,Strain0+2),(D(St:ed,Stress0+1)-D(St:ed,Stress0+2))/1000,'color','b','linewidth',1.5);

xlabel('Axial Strain','Fontname','Times New Roman','Fontsize',16);
ylabel('\sigma_1-\sigma_3 (kPa)','Fontname','Times New Roman','Fontsize',16);
set(gca,'Fontname','Times New Roman','FontSize',16,'linewidth',1.0);

set (gcf,'Position',[100,150,900,400], 'color','w');